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TRACK  INITIALIZATION  SENSITIVITY  IN  CLUTTER 

1.  INTRODUCTION 


Track  initialization  error  is  quantified  as  a  function  of  clutter  density  for  the  nearest  neighbor 
assignment  function  (NNAF).  Quality  of  track  initialization  is  thereby  related  to  assignment 
accuracy.  Two  statistical  distributions  of  total  track  initialization  error  are  derived  under  idealized 
hypotheses  representing  excellent  and  poor  performance  of  the  NNAF  in  clutter.  A  track 
initialization  sensitivity  region  is  obtained  by  constraining  the  mean  initialization  error  under  the 
excellent  NNAF  performance  hypothesis  to  have  (lower  tail)  significance  level  at  most  a  under 
the  poor  NNAF  performance  hypothesis.  The  initialization  error  region  is  derived  for  general 
nonlinear  state  equations  and  for  measurement  equations  with  additive  Gaussian  noise;  the 
initialization  error  region  is  not  conditioned  on  particular  measurement  sets,  but  it  is  a  function  of 
clutter  density  and  the  number  of  scans  in  the  track  initialization.  As  the  number  of  scans  in  the 
initialization  increases,  examples  show  that  sensitivity  to  errors  in  initial  velocity  increases  linearly 
compared  with  sensitivity  to  positional  errors  when  constant  velocity  models  hold. 

There  exists  a  critical  clutter  density  at  which  the  initialization  error  region  is  the  empty  set. 
One  interpretation  of  this  result  is  that  even  perfect  track  initialization  is  a  poor  initialization  in 
clutter  whose  density  exceeds  critical  density.  The  critical  clutter  density,  normalized  by  the 
determinant  of  the  additive  measurement  noise  process,  depends  only  on  the  dimension  of  the 
measurement  vector,  moreover,  it  is  independent  of  the  state  and  measurement  process  models. 
An  explicit  expression  for  the  critical  clutter  density  is  derived. 

The  methodology  is  used  to  define  a  family  of  curves  called  initialization  clutter  envelop 
(ICE)  curves.  ICE  curves  plot  maximum  clutter  density  against  maximum  initialization  error,  and 
they  are  useful  for  estimating  the  maximum  clutter  density  at  which  successful  tracking  can  occur 
for  a  given  initialization  error.  Further  details  are  discussed  in  section  2. 

In  practice,  measurement  assignment  functions  may  use  a  variety  of  application-specific 
“clues”  to  enhance  performance  and  reduce  track  estimation  error.  The  problem  studied  here  is 
restricted  to  assignment  functions  that  are  strictly  kinematical.  Thus,  only  the  initial  target  track 
and  the  target  measurement  process  are  known  to  the  NNAF.  Performance  limits  for  the 
kinematic  problem  are  insightful  and  useful  in  themselves,  and  they  facilitate  quantitative  analysis 
of  the  utility  of  additional  clues  to  the  tracking  process. 
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2.  INITIALIZATION  SENSITIVITY  REGION 


Precisely  one  target  is  assumed  present,  and  an  initial  track  estimate  is  given.  Multiple 
measurements  from  a  sensor  at  T  >  1  successive  independent  scans  are  available  for  re-estimating 
the  target  track.  The  probability  of  detecting  the  target  at  each  scan  is  one,  so  precisely  one 
measurement  comes  from  the  target  in  each  scan.  Measurements  not  originating  from  the  target 
are  termed  clutter.  An  automated  assignment  function  selects  the  measurement  corresponding  to 
the  target  from  among  the  available  measurements  at  each  scan.  The  measurement  assigned  to 
the  target  is  presented  to  the  tracking  algorithm,  and  a  new  estimate  of  target  track  is  computed. 

The  assignment  function  considered  here  is  the  NNAF.  The  NNAF  is  given  a  set  of 
(unlabeled)  measurements  comprising  both  target  and  clutter,  and  the  measurement  closest  to  the 
expected  measurement,  conditioned  on  the  initial  track,  is  assigned  to  the  target  at  each  scan.  The 
NNAF  has  the  opportunity  to  commit  errors,  that  is,  to  select  a  clutter  measurement  as  the  target 
measurement.  An  explicit  mathematical  definition  of  the  NNAF  in  terms  of  measurements  and 
initial  target  track  is  postponed  momentarily  (see  equation  (1)). 

A  null  hypothesis  is  used  to  model  poor  track  initialization.  The  target  is  always  present,  by 
definition,  so  the  null  hypothesis  cannot  be  a  traditional  target  absent  condition.  The  null  must 
pertain  to  the  initial  track  estimate  and  the  strictly  kinematic  information  used  by  the  assignment 
function.  Intuitively,  the  worst  track  initializations  are  so  remote  from  the  target  track  at  every 
scan  that  the  presence  of  any  amount  of  clutter  makes  the  assignment  function  fail  to  select  the 
correct  measurement  with  probability  very  nearly  one.  Idealizing  this  situation  leads  to  the 
following  null  hypothesis  for  poor  track  initialization: 

“7^  :  The  NNAF  never  assigns  the  correct  measurement  to  the  target  track;  that  is,  track 
initialization  is  such  that,  at  each  scan,  the  measurement  assigned  to  target  is  a  clutter 
point. 

Clutter-induced  track  estimation  error  is  the  only  source  of  estimation  error  under  the  null 
hypothesis  Intuitively,  the  likelihood  of  the  NNAF  failing  completely,  as  stipulated  in 
hypothesis  Tf.  diminishes  with  decreasing  clutter  density. 

An  alternative  hypothesis  is  used  to  model  good  track  initialization.  A  simple  model  is  that 
track  initialization  is  so  good  that  the  assignment  function  is  nearly  always  correct  at  any  clutter 
density.  Idealizing  this  situation  leads  to  the  following  alternative  hypothesis: 

X.  The  NNAF  is  error  free;  that  is,  track  initialization  is  such  that,  at  each  scan,  the 
measurement  assigned  to  the  target  is  the  correct  measurement. 

Target  measurement  error  is  the  only  source  of  track  estimation  error  under  hypothesis  "K. 
Intuitively,  the  likelihood  that  the  NNAF  never  makes  a  wrong  assignment,  as  required  by 
hypothesis  X,  diminishes  with  increasing  clutter  density. 
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The  statistic  of  interest  under  1?^  and  :?ris  a  normalized  innovation  statistic,  denoted  S,  which 
measures  the  total  squared  error  (relative  to  the  measurement  covariance)  of  the  given  track 
initialization.  In  sufficiently  heavy  clutter,  the  assignment  functions  modeled  by  and  ^  are 
difficult  to  distinguish  using  the  statistic  S.  When  clutter  reaches  a  certain  critical  density,  the 
nearest  neighbor  measurement  and  the  target  measurement  are,  on  average,  equally  close  to  the 
expected  measurement  (conditioned  on  the  initial  track).  Consequently,  in  clutter  of  density 
greater  than  or  equal  to  critical  density,  the  target  track  is  lost  against  'the  noisy  background 
despite  the  assumption  that  a  target  is  present.  Equating  the  means  of  the  probability  density 
functions  (PDFs)  of  S  under  ^  and  ^  gives  a  fundamental  equation  relating  the  clutter  density  to 
track  initialization  and  provides  a  convenient  definition  of  a  critical  clutter  density  (see  section  3). 

Standard  approximations  to  the  PDFs  of  S  under  and  "K  are  derived  and  used  in  this  report. 
These  approximations  are  exact  under  the  following  modified  hypotheses; 

The  NNAF  is  given  only  clutter  data,  i.e.,  target  is  absent. 

The  NNAF  is  given  only  target  data,  i.e.,  clutter  is  absent. 

The  modified  hypotheses  are  appropriate  for  detecting  target  absence  or  presence.  They  are 
adequate  approximations  for  the  purposes  of  studying  track  initialization  sensitivity.  Moreover, 
they  provide  significant  insight  into  the  problem.  For  these  reasons,  the  hypotheses  '  and  ' 
are  studied  in  this  report.  Numerical  investigation  of  track  initialization  under  hypotheses  and  TC 
will  be  the  subject  of  a  later  study. 

Distributions  for  the  statistic  S  under  '  and  ^ '  differ  from  those  under  1?  and  X,  a  fact  that 
went  completely  unnoticed  in  the  literature  until  pointed  out  in  a  recent  paper  by  Li.*  Under 
hypotheses  and  K,  the  input  to  the  NNAF  comprises  unlabeled  target  and  clutter  measurements, 
yet  the  NNAF  does  not  happen  to  select  incorrect  measurements,  namely,  the  target  measurement 
under  ^  and  a  clutter  measurement  under  >?r.  Given  the  opportunity  to  commit  error,  hypotheses 
and  state  that  no  errors  occurred.  In  contrast,  under  hypotheses  '  and  X ',  the  NNAF 
cannot  select  incorrect  measurements  because  the  inputs  to  the  NNAF  are  censored;  the  target 
measurement  is  absent  under  '  and  clutter  measurements  are  absent  under  ^ '. 

A  track  error  sensitivity  region,  denoted  by  ^'(a) ,  is  defined  by  setting  a  significance  level  a 
for  rejecting  '  and  requiring  the  mean  of  the  statistic  S  under  ^ '  to  be  significant  at  this  level, 
i.e.,  the  mean  of  S  under  ^ '  lies  on  the  lower  tail  of  S  under  1?'' '  at  the  point  of  significance  a. 

The  defining  inequality  of  the  initialization  region  is  derived  in  this  section,  as  well  as 

receiver  operating  characteristic  (ROC)  curves  under  hypotheses'^ '  and 

A  sensor  measurement  is  a  vector  of  fn>\  real-valued  components;  that  is,  each  sensor 
measurement  is  a  point  in  R"" .  The  normalized  distance  between  two  measurements,  say  z  and  w, 
is  defined  by 
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d{z,w)  =  , 


(1) 


where,  throughout  this  report,  R  is  the  covariance  matrix  of  additive  Gaussian  measurement  noise 
(see  equation  (8)).  The  normalized  distance  (equation  (1))  is  used  by  the  NNAF.  The  sensor’s 
field  of  view  is  assumed  to  comprise  all  of  ,  a  simplifying  idealization.  Clutter  is  modeled  as 

uniform  Poisson  (point)  clutter  over  /?”  with  normalized  density  parameter  ,  w^here  \R\  is 
the  determinant  of  R.  The  normalized  density  is  dimensionless;  this  is  a  practical  advantage  for 
application  comparisons.  A  theoretical  justification  for  the  normalized  density  is  as  follows:  let 
R  =  LL' ,  where  L  is  the  Cholesky  factor  of  R.  Then  the  transformation  z  =  Ly  transforms  the 

normalized  distance  (equation  (1))  into  the  usual  Euclidean  distance  in  /?'”  and,  simultaneously, 

I 

transforms  the  Poisson  clutter  process  into  another  Poisson  process  with  density  A|  |  Z,“’ j  =  . 
That  the  density  parameter  transforms  in  this  way  is  an  immediate  consequence  of  the  Mapping 
Theorem  for  Poisson  processes  (see  Kingman,  section  2.3).^ 

Let  denote  a  normalized  sphere  in  R"  with  radius  r  and  centered  at  any  fixed  point, 

say  the  origin,  that  is, 

0^{r)  =  [z  GR'”:d{z,0)<r} . 

The  volume  of  (r)  is^ 


IRl^r" 

r(f+i) 


-  \W  r" 


(2) 


Using  properties  of  the  Gamma  function  and  r\—\^n^  gives  the  special  cases 


r.  An  ,  n 

K.  =2,  n,  K.  ~  — ,  and  k.  =  — . 
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It  is  well  knowm"*  that  the  PDF  of  the  Euclidean  distance  of  the  Ath  closest  clutter  point  to  the 
origin  (or  any  fixed  point)  is 


r 


m\ 


l^ 


P,(r)  = 


V _ J 

(k-l)! 


r*”  ’  exp\ 


^  i  ^ 

-KAWr 
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Only  the  PDF  of  the  closest  clutter  point  is  needed  in  this  application,  and  this  is  the  case  k  =  \  ', 
explicitly, 

pM  =  exp{-K„X\R\^r'”). 

The  PDF  of  the  squared  distance  of  the  closest  clutter  point  is,  therefore, 

where  u  =  r^ .  The  moments  of p(u)  are 


2v 


-  ^CO 

u'’  =  u''p(u)  du  = 


m+i)  _-r(¥+i) 


7t 


■ra+i) 

X\R^ 


V  =  0,1,2,..., 


1  m 


as  is  seen  via  the  change  of  variables  ^  ^  ^ .  In  particular. 


TT 


X\R^ 


(4) 


is  the  mean-squared  distance  of  the  closest  clutter  point  to  the  origin. 

A  total  of  T  >  1  successive  independent  measurement  scans  are  given.  Let  r,  denote  the 
normalized  distance  (equation  (1))  of  the  closest  clutter  measurement  to  the  origin  at  time  t. 

Under  the  null  hypothesis^ ',  the  squared  distances  m,  =  rf  are  independent,  so  it  follows  that  the 
total  innovation  statistic 


(5) 


has  mean  and  variance 


=  Tu  and  var{S^^  = 


(6) 


The  statistic  S^.  is  independent  of  the  target  initialization  because  Poisson  clutter  is  spatially 
homogeneous  in  the  measurement  space,  R"' .  The  PDF  of  the  statistic  (equation  (5))  will  be 


denoted  by  p^s 


K 


X 

^ ,  T,  ‘J5*’') .  No  explicit  form  for  this  PDF  is  available  except  for 


m 


=  2; 
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however,  the  appendix  shows  that  useful  thresholds  for  hypothesis  testing  purposes  are  readily 
obtained  for  all  m. 

Let  R"  denote  the  target  state  space.  A  target  initialization  comprises  a  sequence  of  target 
states  of  lengthT,  denoted  where  each  state  gR”.  Position¬ 

keeping  (PK)  track  initialization  is  assumed;  i.e., 

t  =  (7) 

PK  initialization  is  thus  fully  characterized  by  the  first  point  .  Additive  Gaussian  noise 
models  are  assumed  for  the  given  nonlinear  target  measurement  process;  hence,  the  target 
measurement  process  is  modeled  by 

z,^N,(xr)  +  «’,,  (8) 

where  the  additive  noises  }  are  independent  and  identically  distributed  with  covariance  matrix 

R  and  denotes  the  actual  (but  unknown)  target  states.  Under  hypothesis  X ',  the  correct 

target  measurement  is  known — only  the  initialization  is  erroneous.  Hence,  the  measurement  error 
conditioned  on  the  initialization  is  given  by 

= z,  -  H,{xr) + HXxT)  - 
=w,+HXx'r)-H,ixr^). 

The  errors  {£■,  }  are  Gaussian  distributed  with  mean  vector  and  covariance 

matrix  R.  The  statistic 

s,,  =  i;(z,-//,(xr"))  (*;■“))  (9) 

t=\ 

is  distributed  =  Zmii'  1^) ,  that  is,  noncentral  chi-squared  with  mT  degrees  of  freedom  and 

noncentrality  parameter 

•5 = staL”')-  -ff.l*:")) «-'(/?, L'”*)-  //.(x")).  (10) 

/=:! 

This  result  follows  easily  from  Muirhead  (Theorems  1 .3.4  and  1.4. 1).^  The  noncentrality 
parameter  is  the  total  squared  initialization  error.  The  mean  and  variance  of  ,  are  (see 
Muirhead,  page  24f 


7 


(11) 


S^,  =mT+6  and  var{^S^.  )-2mT+4S. 


Because  ofPK  initialization,  the  statistic  S^,  is  conditioned  solely  on  the  first  point  x 

^initial 


initial  r» 

1  of 


If  measurement  process  is  linear,  so  that  =  for  all  x,  then  the  total  squared 

initialization  error  can  be  written  in  the  form 


t=\ 

Further,  if  the  target  process  is  linear,  so  that  =  F,x  for  all  x,  then 
d=A'  WA, 


where 

A  (13) 

is  the  track  initialization  error  vector  at  time  t=  I,  and  the  matrix  W\s  given  by 


T  (  ^  f  t-\ 

f=l  V  >=0  J  V  ;=0  J 


where  is  the  n  x  77  identity  matrix.  If  the  target  process  is  stationary  also,  so  that  ^^(x)  -Fx 
for  all  t,  then  W  simplifies  further  to 


f‘~')  i?-'  (f,  f‘-'). 


(14) 


Further  simplification  of  W  seems  possible  only  in  special  cases.  The  matrix  fV  is  clearly  positive 
semidefinite,  as  is  evident  from  equation  (12).  In  many  practical  problems,  F  and  }  are  such 
that  Wis  positive  definite;  however,  no  general  condition  for  the  positive  definiteness  of  IF  is 
presented  here. 

Let  the  significance  level  a  for  rejecting  the  hypothesis  ^ '  be  specified.  If  the  normalized 

clutter  density  is  such  that  ,  where  is  defined  in  section  3,  then  the 

quantity  s^{k „X\P\^ ,T)  satisfying  the  equation 


8 


1-a 


!  Ps(^ 

i 


(15) 


defines  the  threshold  at  significance  level  a  for  rejecting  '  in  a  traditional  lower-tail  significance 

1 

test.  The  threshold  s^{k ,T)  depends  on  a,  the  normalized  clutter  density,  and  the  number 
of  scans  T  contributing  to  the  innovation  statistic.  It  follows  by  the  change  of  variables 


^=[K,x\4]'n 


that 


-2 


-  f  - 

V  J 


where  the  standardized  threshold  is  defined  by 


(16) 


1-a  =  J  P^(7J  jl,  r, .  (17) 

Computation  of  the  standardized  threshold  ■S'q.(T)  is  discussed  in  the  appendix,  where  several 
tables  are  also  given. 

The  initialization  sensitivity  region  of  significance  a,  denoted  ,  is  defined  by 

constraining  the  mean  of  the  innovation  statistic  under  '  to  be  less  than  or  equal  to  the  threshold 

y 

s^{k ,T) .  The  mean  is  the  noncentrality  parameter  (equation  (10)),  so  the  region  ;?'(«)  is 
given  by  the  quadratic  inequality 


(18) 

where  the  identity  (equation  (16))  is  substituted  in  equation  (18).  If  IT  is  positive  definite,  the 
significance  region  is  a  nondegenerate  ellipsoid  in  the  target  state  space,  R" .  If  IT  is  not 

positive  definite,  the  significance  region  has  infinite  extent  in  directions  corresponding  to  the 
eigenvectors  of  the  zero  eigenvalues  of  IT.  Track  initialization  in  the  significance  region  ^'[a) 

guarantees  that  the  mean  of  S^.  does  not  exceed  the  significance  level  a,  assuming  PK  track 
initialization.  Intuitively,  initializing  in  ^'(a)  guarantees  that  at  most  a  percent  of  all  potential 
tracks  generated  by  the  NNAF  in  clutter  are  “track-like”  in  comparison  with  real  tracks. 


^’{a) 


zl  eR"-.A’WA  < 


-2 
Am 


'.(T) 
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'  A  more  conservative  initialization  region  is  defined  by  specifying  an  additional  significance 
level  on  the  alternative  hypothesis  "K ' .  Let  /?  denote  the  lower-tail  significance  level  for  the 
statistic  ,  that  is,  1  -  /?  is  the  probability  of  rejecting  TC.' .  The  initialization  sensitivity  region  at 
the  operating  point  {a,P)^  denoted  1R'{a,0),  is  defined  by 


Jei?" 


J 


A'  WA) 


(19) 


For  every  probability  a  there  exists  a  probability  such  that 
;R'[a,0)Q^'{a),foT  ally^<y^o, 

that  is,  the  initialization  region  is  smaller  than  the  region  ^'(a)  for  all  /?  <  .  A 

rough  approximation  is  ^'(a) «  • 


The  ROC  curve  for  '  and  ^ '  is  the  locus  of  the  point 


(20) 


as  s  ranges  from  zero  to  infinity.  The  integrals  range  over  the  interval  (0,  5)  instead  of  the  more 
typical  interval  (•s.oo)  because  both  hypothesis  tests  are  lower-tail  tests.  The  ROC  curves 
(equation  (20))  are  functions  of  clutter  density,  initialization  error,  and  T,  so  they  are  actually  non- 
planar  surfaces  and,  hence,  difficult  to  visualize. 

Insightful  planar  curves  are  obtained  in  the  following  manner:  set  the  operating  point  {pc,P) 
and,  for  every  threshold  s,  define 


S\s)  =  max\5\P  <  ^ \ 


(21) 


X\^)\eA  =max\x\BA  -.a  <J  |  X\RX  ,T,'»')d^ 


(22) 


It  can  be  verified  that  the  equation  defining  the  set  in  equation  (21)  has  a  unique  solution  for  each 
5  >  s(p) ,  where  s(P)  is  the  p  percent  point  of  the  (central)  chi-square  distribution,  i.e.. 
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(23) 


•in) 

yff=  I  I  m- 

0 

For  s<s{0),  the  set  in  equation  (21)  is  empty  because  no  solution  to  the  defining  equation  exists. 

The  curve  (.^^’(5),/T(5)|i?|2)  defined  parametrically  as  s  ranges  from^O^  to  infinity  is  called  the 

ICE  curve,  and  it  depends  only  on  the  number  of  scans  T.  ICE  curves  are  useful  in  approxi¬ 
mating  the  maximum  density  clutter  in  which  targets  can  be  successfully  tracked.  Knowledge  of 
the  sensor  signal  processor  and  clutter  density  determines  target  signal-to-noise  ratio  (SNR),  so 
maximum  clutter  density  estimates  are  equivalent  to  minimum  target  SNR  estimates. 


3.  CRITICAL  CLUTTER  DENSITY 


As  discussed  in  section  2,  critical  clutter  density  occurs  when  the  first  moments  of  the 
statistics  and  are  equal.  Using  the  hypotheses  ’  and  ^ '  and  equating  the  mean  values  of 

the  statistics  S^,  and  S^.  gives 


S^,  =  T  u  =  mT + S  =  S^, 


(24) 


Solving  for  the  total  squared  initialization  error  S  gives  the  fundamental  equation  relating  clutter 
density  and  track  initialization  error,  namely. 


S  =  t{u  -m^=  rp(/l.)  >  0, 


(25) 


where 


p(X)  = 


r  1 

2 

^(i+1) 

/’(f+i) 

m 

-  m 

TV 

_  A|/?P  _ 

(26) 


The  clutter  density  parameter  A  appears  only  on  the  right-hand  side  of  the  equation.  Solving  the 
inequality  p{^)  ^0  for  A  gives 


A< A  ..  ,  = 

~  critical 


r(i+i) 

Tvm 


m 

~2 


(27) 
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The  critical  clutter  density  is  a  function  only  of  the  measurement  dimension  and  the  determinant  of 
the  covariance  matrix  of  the  additive  Gaussian  noise.  The  normalized  critical  clutter  density, 

1 

^critica\^'^  depends  only  on  the  measurement  dimension,  a  fact  that  is  intuitively  satisfying. 

The  critical  clutter  density  is  independent  of  the  number  of  scans  T,  as  might  be  anticipated 
from  the  problem  formulation.  More  importantly,  it  also  holds  for  general  nonlinear  target  and 
measurement  models.  Different  target  and  measurement  process  models  (i.e.,  the  functions 
and  {^,} )  yield  the  same  critical  clutter  density.  This  is  not  to  say,  however,  that  different 
models  track  equally  well  in  the  same  density  clutter. 

The  critical  clutter  density  is  defined  via  an  extreme  circumstance,  namely  when  the  means  of 
two  PDFs  are  equal.  As  such,  the  critical  density  should  be  interpreted  as  an  upper  bound  at 
which  a  strictly  kinematic  tracking  algorithm  can  be  expected  to  operate.  Such  a  performance 
limit  provides  a  quantitative  standard  of  comparison  for  alternative  tracking  algorithms  in  clutter. 
For  example,  a  tracker  that  functions  reliably  at  normalized  clutter  densities  up  to  10  percent  of 

^ critical  ^sensor  iP  probabty  better  than  one  that  functions  only  up  to  1  percent  of 

j_ 

^criticai\Kemor  2^  '  ^ho  Utility  of  the  normalized  clutter  density  as  a  standard  comparison  measure 
deserves  further  examination  in  practical  problems. 

Particular  values  of  the  normalized  critical  density  are,  using  equation  (27), 


X 


critical 


W  = 


S’ 


_l_ 

Itv' 


A 

18;r 


HI 


1 

327r 


for  measurement  vector  dimensions  m=  1,2,  3,  and  4,  respectively.  Using  Stirling’s 
approximation  to  the  Gamma  function  and  the  curious  limit 


(28) 


+  =  e  ^  =0.561459, 


(29) 


where  y  =  0.577216  is  Euler’s  constant,  gives  the  asymptotic  result 


^ critical  “ 


^  e  '  -JmTV 


{27isY 


as  w  ->  00. 


(30) 
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It  follows  from  equation  (30)  that 


Urn 


critical  m 


critical  m+1 


=  .I^  =  A.\3TJ. 


(31) 


Evidently,  the  critical  clutter  density  decreases  with  increasing  measurement  dimension  m  by  a 
factor  of  about  four  for  additional  components  of  the  measurement  vector.  As  shown  by  the 
tabular  data  below,  these  limiting  values  are  approached  rather  slowly. 


m 


1 

0.70711 

4.4429 

2.9365 

2 

0.15916 

4.0387 

1.9315 

3 

0.39407  (-1) 

3.9617 

1.6137 

4 

0.99472  (-2) 

3.9472 

1.4549 

5 

0.25201  1-2) 

3.9504 

1.3653 

10 

0.25589  (-5) 

3.9965 

1.1817 

15 

0.24660  (-8) 

4.0287 

1.1210 

20 

0.22981  (-11) 

4.0492 

1.0907 

25 

0.20958  (-14) 

4.0631 

1.0725 

100 

0.24098  (-60) 

4.1130 

1.0181 

1000 

0.18260  (-614) 

4.1307 

1.0018 
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4.  EXAMPLE:  CONSTANT  VELOCITY  TARGET 


The  example  is  uses  a  constant  velocity  target  model  with  two  dimensional  measurements. 
For  ease  of  exposition,  measurement  vector  components  are  referred  to  as  position  coordinates. 
Measurement  scans  are  available  at  fixed  sampling  intervals  of  r  >  0.  The  units  of  time  and 
position  (e  g.,  seconds  and  meters)  are  irrelevant  because  standardized  error  coordinates 
(described  below)  are  used  throughout.  Extension  to  constant  acceleration  target  models  is 
straightforward  and  not  pursued  in  this  report. 


In  the  example  m-2,n  =  4,  and  the  constant  velocity  motion  model  is 

1  T 

0  1 
F  = 

0  0 
0  0 

The  target  state  components  are  position  x  and  velocity  x .  From  the  definition  (equation  (14)), 
IF  is  a  4x4  block  diagonal  matrix  with  the  2x2  diagonal  blocks 


0  0 
0  0 
1  T 
0  1 


H  = 


10  0  0 
0  0  10 


and  R  = 


<7 

0 


(32) 


T 

\tT{T-\) 


\tT{T-\) 

\FT{T-\){2T-\) 


(33) 


is  singular  for  f  =  0  because,  in  this  case,  the  matrix  F  models  a  motionless  target.  For 
T>0,  the  matrix  is  singular  for  T  =  1  because  velocity  is  unobservable  from  a  single  position 
measurement.  For  the  remainder  of  this  section,  it  is  assumed  that  r  >  0  and  T>2 .  Let 
(x,x,  >/,>')  denote  the  (transposed)  target  state  vector.  Standardized  error  coordinates  for 
position  and  velocity  are 

f  Ax  tAx  Ay  r4y^ 

V  O'  ’  o’  o’  o  / 


1  /A  A‘  A  A*\  /,..true  ^initial  .^true  •  initial  ..true  initial  ,*,true  •  initial  \  'T'U... 

where  (^,  4br,  4y,  4y)  =  (x,  -x,  ,  x,  -x,  .  y,  -y,  ,  y,  -y,  ).  The  block 

structure  of  W  and  the  explicit  form  (equation  (33))  gives  the  initialization  region 


/ 

^Ax^ 

/ 

a 

Ax 

w° 

a 

Ax 

+ 

0 

Ay 

w° 

0 

Ay 

V  cr^ 

^  a  y 

where  the  standardized  matrix  W°  is  defined  by 


(34) 
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jr(r-i) 


Because  the  x  andj  coordinates  separate  in  (34),  the  initialization  region  in  this  case  is  most  easily 
studied  by  supposing  no  initialization  error  in  one  component.  The  eigenvalues  of  W°  are 


For  r  >  2 ,  it  is  easily  seen  from  these  forms  that  the  eigenvalues  are  distinct  and  positive. 
Consequently,  ,  Wd  ,and  W  are  positive  definite. 

The  positive  branch  of  the  square  root  in  equation  (35)  gives  the  larger  of  the  two 
eigenvalues.  The  lengths  of  the  major  and  minor  axes  of  the  elliptical  track  initialization  region 
^ '  in  standardized  error  coordinates,  are  proportional  to  the  square  roots  of  the  smallest  and 

largest  eigenvalues,  respectively,  of  W° .  Thus,  the  eccentricity  of  the  ellipse  is 


eccentricity  = 


major  axis  length 
minor  axis  length 


For  7’=  2,  5,  10,  15,  20,  and  25,  the  eccentricity  rounded  to  the  nearest  integer  is  3,  5,  10,  16,  22, 
and  27,  respectively.  Asymptotically, 

+o(r')  and  y_=\T+o(\),  asr->oo,  (36) 

so  that 

eccentricity  ^  T  +  o(l),  as  T  — »  oo.  (37) 

The  initialization  region  eccentricity  therefore  increases  linearly  with  T  for  large  T. 

The  eigenvectors  of  W°  corresponding  to  are  given  by 


16 


where 


<>.= 


T-l 


-l  +  K7'-lX2r-l)±Vl  +  i(7’-lX7’-2)+i(?'-l)'(27--l)= 


(39) 


The  angle  between  the  positive  standardized  position  axis  and  the  eigenvector  e_  corresponding 
to  major  axis  of  W°  is  =  tan  .  For  7=2,  5,  10,  15,  20,  and  25,  y/  is  -58°,  -19°,  -9°,  -6°, 
-4°,  and  -3.5°,  respectively.  It  is  easily  verified  that  6  <  0  for  7  >  2,  so  y/  must  be  negative. 
Asymptotically, 


For  sufficiently  large  7,  the  angular  “tilt”  of  the  initialization  region  diminishes  linearly  with 
increasing  7.  The  negative  correlation  between  normalized  position  and  velocity  errors  accords 
well  with  intuition. 

The  half-lengths  and  L_  of  the  major  and  minor  axes,  respectively,  of  the  region  are 
functions  of  clutter  density  and  significance  level  a.  Explicitly, 


(41) 


where  the  standardized  threshold  is  tabulated  in  the  appendix.  From  equation  (28),  the 

critical  clutter  density  for  two-dimensional  measurements  is  — ,  so  a  normalized  clutter  density 

In 

representative  of  moderate  to  heavy  clutter  is 


i 


2n 


For  7=  2,  5,  10,  15,  20,  and  25  and  for  a  1-percent  significance  level  (i.e.,  or  -  0.01),  the  major 
axis  half-lengths  are 

=  2.789,  4.140,  5.409,  6.044,  6.434,  and  6.704, 
respectively;  the  minor  axis  half-lengths  are 

L_  =  1.065,  0.8737,  0.5317,  0.3818,  0.2986,  and  0.2457, 
respectively.  These  half-lengths  assume  no  error  in  one  component.  For  errors  in  both  x  and>' 
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components,  a  more  appropriate  “average”  initialization  region  is  obtained  by  dividing  these  half- 
lengths  by  Vs  ,  assuming  equal  error  contributions  in  each  component. 

Figure  1  shows  the  initialization  region  ^'(a)  corresponding  to  1 -percent  significance  levels, 
respectively,  for  rejecting  clutter  for  several  choices  of  batch  length  T.  The  relative  sizes  of  the 
major  and  minor  axes  of  ^’(a) ,  as  well  as  the  orientation  angle  y/ ,  agree  with  the  results  given 

above. 

It  is  evident  from  figure  1  that,  as  T  increases,  track  initialization  error  in  standardized 
coordinates  is  more  sensitive  to  initial  velocity  errors  than  to  initial  position  errors.  This  trend 
confirms  the  asymptotic  formulae  given  above.  It  follows  that  the  allowable  error  in  Ax  (and  in 
Ay )  is  inversely  proportional  to  the  sample  interval  r. 

ICE  curves  are  described  in  section  2  as  a  way  to  quantify  the  following  intuitively  reasonable 
statement:  increasing  the  normalized  clutter  density  requires  progressively  better  track 
initialization  to  maintain  a  fixed  significance  level  for  rejecting  clutter.  However,  ICE  curves  are 
not  presented  here. 


tAx 


Figure  1.  Initialization  Error  Regions  for  Significance  Level  a  =  0.01 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


Further  study  of  several  topics  is  recommended.  The  utility  of  families  of  ICE  curves  as  a 
means  of  quantifying  the  maximum  clutter  density  in  which  target  tracking  can  be  reliably 
performed  merits  further  study.  The  development  of  initialization  error  regions  for  constant 
acceleration  target  models  also  merits  further  study.  The  correct  hypotheses  ^  and  X  remain  to 
be  investigated.  Finally,  when  non-kinematic  information  is  exploited  by  the  assignment  function, 
the  approach  of  this  report  is  inadequate  and  must  be  reformulated.  The  best  way  to  approach 
this  topic  depends  on  the  application  and  the  nature  of  the  non-kinematic  data  used. 

The  approach  developed  is  applicable  to  any  filtering  problem  in  which  initialization  in  clutter 
is  potentially  a  problem.  The  example  presented  suggests  that  this  formulation  provides  a 
reasonable  and  useful  method  to  assess  the  relative  sensitivity  of  different  state  variables  to 
initialization  in  clutter. 
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APPENDIX 

COMPUTATION  OF  SIGNIFICANCE  THRESHOLDS  UNDER 

The  standardized  threshold  5^(7)  is  defined  by  equations  (16)  and  (17).  For  m  =  2,\X  is 
easily  shown  that 

p\rj[\,  T,  = 

so  that  the  cumulative  function  is,  in  this  case,  the  incomplete  Gamma  function  (see  Abramowitz 
and  Stegun,  equation  (6.5.1)).®  Algorithms  for  computing  the  incomplete  Gamma  function  are 
widely  available  (e.g.,  in  Mathematica),  so  computing  the  threshold  is  simply  a  matter  of  solving 
the  nonlinear  equation  (17).  \n  Mathematica,  Newton-Raphson’s  iteration  for  the  solution  is  the 
following  one  line  expression; 

FindRoot  [GammaRegularized[T,  s]= 1  -alpha,  {  s,  sO }  ] ,  ( A-2) 

where  50  is  an  initial  guess  at  the  solution,  and  GammaRegularized{  ]  is  an  intrinsic  Mathematica 
function  that  equals  one  minus  the  incomplete  Gamma  function.  The  algorithm  (equation  (A-2)) 
is  insensitive  to  the  choice  of  initial  point  50. 

For  w  ^  2 ,  the  required  PDF  is  not  available  in  closed  form,  so  computing  the  thresholds 
5„(7’)  accurately  requires  numerically  evaluating  the  characteristic  function  f{^)  oipiu)  and 

integrating  \f(^^ .  For  full  details,  see  Nuttall.’  For  the  present  purposes,  however, 
approximate  thresholds  are  easily  computed  by  simulation  using  the  following  procedure.  The 

exceedance  of p{u)  for  =  1  is,  from  equation  (3), 


E(v)  -  '  exp[  -u ^  JJm  =  exp 


-V 


(A-3) 


Its  inverse  function  is 


2 

V  =  (-/«£)  . 


(A-4) 


Hence,  ifE  is  uniformly  distributed  on  the  interval  (0,1),  then  v  is  distributed  according  Xop{u). 
Consequently,  the  sum 


T 


A-1 


is  distributed  |l  ,  provided  each  term  v,  is  computed  via  equation  (A-4).  Samples 
of  the  sum  (equation  (A-5))  are  generated  using  equation  (A-4)  and  a  pseudo-random 
number  generator  and  sorted  so  that  -^(i)  -  -^(2)  where  denotes  the  «th-order 

statistic.  The  lower-tail  significance  point  s^(T)  is  estimated  as 

where  [oiV]  in  equation  (A-6)  denotes  the  greatest  integer  less  than  or  equal  to  ccN .  Sample  size 

N  strongly  affects  the  accuracy  of  the  estimates  (equation  (A-6)),  especially  for  small  a .  For 
sufficiently  small  values  of  a,  the  methods  of  reference  7  are  recommended  over  simulation. 

The  tables  below  for  m^2  are  compiled  using  #=  10,000  samples.  These  thresholds  are 
probably  accurate  to  about  two  digits.  For  m  =  2,  the  numerical  procedure  described  above  is 
used  to  obtain  thresholds  that  are  rounded  to  the  precision  shown  in  tables  A-1  and  A-2. 


Table  A-L  Standardized  Thresholds  ^^{T)  for  m  =  1  and  m-  2. 


A-2 


Table  A-2.  Standardized  Thresholds  (?’)  for  m-  3  and  m-4 
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